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A generalization of Kermack-McKendick model of epidemics to the case of inhomogeneous susceptibil- 
ity of population is proposed. Some quantitative and qualitative features of epidemic process development 
in this situation are established. 

The propagation of epidemics in a biological population (collective) is a complex pro- 
cess. In a widely used Kermak-McKendric model of epidemics [1] (also known as SIR- 
model), and in many others, the basic characteristic of process, the susceptibility coeffi- 
cient A, is supposed to be constant over a population, i.e. all individuals are supposed 
to have the same immunity. Numerous researches, however, show that this assumption 
is not satisfied. For various reasons wide fluctuations of susceptibility parameters are ob- 
served in real populations, and they may influence significantly the dynamics of epidemic 
process. In the present work we generalize the SIR-model to the case of inhomogeneous 
susceptibility (immunity) of the population. 

To obtain the equations of the model we shall assume that a closed population of N 
individuals can be split into n groups (2 < n < N), such that all individuals from the i-th. 
group have the same susceptibility coefficient Aj > 0, i = l,n. The system of equations of 
Kermack-McKendrick for this case has the form (the dot stands for the time derivative): 

n 

5, = -W, 7 = 7^A^- 7 7, R = -yI, (1) 

i=l 

where Si = Si(t) is the number of individuals in the i-th group susceptible to an infection 
at the moment t, I — I(t) is the number of infected at the moment t, R — R(t) is the 
number of recovered at the moment t, and 7 > is the " factor of elimination". The 
initial conditions for the system (1) have the form: 

5,(0) = S 0t , J(O) = Jo, R(0 = 0. (2) 

Notice that in model (1), similar to elementary models of the chemical kinetics, we 
neglect the spatial distribution of individuals. 

To estimate the effect of inhomogeneity of the immunity we are going to compare the 
characteristics of the inhomogeneous population under consideration and of the corre- 
sponding homogeneous population with a susceptibility given by 

Ahom — A — , {O) 

•JO 

where So = Ym=i $oi is the total number of susceptible individuals at t = 0. The average 
susceptibility is constant for homogeneous populations, whereas in the inhomogeneous 
case it is easily seen to be a decreasing function of time: 
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which significantly influences the process. In the initial moment we have 

n 

Snonhom(fy ~ S(0) = ~Iq XiS 0i + \IqSq = 0, 

i=l 

n n 
5'?ion/iom(0) — 5*0 = -M^^ A 2 Sio — AjSoi) 2 ] < 0, 

1=1 1=1 



(4) 



(5) 



where S'(O), >S(0), 5' nori / lom (0), S' no „/ lom (0) are the "speed" and "acceleration" of the process 
at t = for the homogeneous and inhomogeneous populations, respectively. 

Thus, at t = the speeds are equal, but the accelerations are different, and \S non h om (0)\ < 
^(O)]. Hence, the inequality (4) implies that at the initial stage the epidemic process in 
the inhomogeneous population develops slower than in the homogeneous one. 

From (1) we have the following "conservation law": 



^e" Al A I(r)dT + I(t) + 7 / I(r)dr = N. 
i=i Jo 



(6) 



Let s* = / °° I(t)cIt. This quantity is related to the size of epidemics: z(oo) = 7s* 
and, under the assumption that I (00) = 0, is a root of the equation 



n 



i=i 

The quantities s* and z(oo) are continuous functions of parameters Aj, and s* varies from 
/0/7 to iV/7 monotonically in Aj. 

We are now going to find approximative (asymptotic) solutions to the transcendent 
equation (7). To this end, we shall consider three ranges of values of parameters. 

1. In the range 

A,s* > 1 (8) 

we have: 

N 

s* — and z(oo) = N, (9) 

7 

that is, the epidemic will finally get the whole population. 

2. In the range 

X iS * ^ 1 (10) 

we get 

s *^I(jV-^°) (11) 

7 e 



and, correspondingly, 



Ai - Ai - . . . - A n ~ ^ So ' ( 12 ) 

e 

2(00) =iV-—. (13) 
e 



3. In the range 



A,s* < 1, (14) 
developing the exponentials in (7) up to the second order terms, we find: 

Hence, in this range, the inhomogeneity of the population can affect the size of epidemics. 
In turn, the size of epidemics in the corresponding homogeneous population with the 
susceptibility A has the form: 

Zhomioo) = = °^ _ . (16) 

7-5 A + v /(7-^oA) 2 + 2/ 5 A 2 

It is not difficult to show using the Silvester criterion that 

z(oo) < z hom (oo), (17) 

that is, the size of epidemics in the homogeneous population is greater than in the inho- 
mogeneous one. 

Let us suppose now that the population is large and the quantity A is a random 
variable. We shall assume that in the limit n — > iV the actual susceptibility of the 
population can be described by some non-negative continuous probability density /(A), 



/(A)dA = l. (18) 







Let 



and 



POO 

E\ = \* = / Xf(X)dX (19) 
Jo 



a 2 = EX 2 - (EX) 2 (20) 

be, respectively, the expectation and the dispersion of susceptibility. The estimate (15) 
then becomes: 

2I °'< (21) 



b(X') + V(MA*)) 2 + 2I a(A',<7 2 )' 
where 
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a(A*) = S (\* 2 + a 2 ), b(\*)= 1 -S \*. (22) 

Notice that the integration in (14) - (15) is taken over the positive real axis, for the case 
A < corresponds to individuals absolutely resistant to infection, which, as a rule, are 
non-existent in real populations. 

Let us find now the influence of the dispersion of susceptibility on the size of epidemics. 
Let A* be fixed. It then follows from expressions (21)-(22), that the size of epidemics grows 
as a 2 decreases and reaches the maximum 

lim z(oo) = °^ , (23) 

^ 7 - SoA* + ^(7 - SoA*) 2 + 2/qSoA* 2 

which is consistent with the discrete case (11). 

Let's pass to the managing parameter p = 7/ A* and we shall analyse the expression 
(23) at Io <C (p — So) 2 . The relation for the size of epidemie then becomes more simple: 

,(00) - -P{P-So -\So-P\] , (24) 

In the considered area of parameter's values two cases are possible: 

1. p < So, that is carried out at enough small Jo, i.e. at small number of ills. In this 
situation 



and 



z(oc) - 2p(l - (25) 
lim z(oo) = 2p. (26) 

So— >oo 

Thus, the size of epidemic in case of the big population with small initial number of ills 
will be limited and to be determined by the quantity of the managing parameter. 
2. So < p. Then for the size of epidemic we shall obtain 

z(oo) = 0. (27) 

So, within the framework of the considered model it was possible to show, that het- 
erogeneity of a population on a susceptibility (immunity) influences on epidemic process. 

Restrictions on the received by us in the given work areas it is possible to interpret, if 
to take into account the expression for average relative speed epidemic "waves": 

Vi = ~ J ° T td \ (28) 

where Tj is the duration of epidemie in i-th group. Then the expression (10), determining 
the second area, it is possible to replace on the expression 

v = ^ (29) 

where T is the duration of epidemie in a population, v is average speed of the epidemie 
"waves". 
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For relative speed, using the first equation of system (1), we shall obtain: 



where 



Then for the first area 



v + A / = 0, (30) 



Ao = ^^, « = (31) 



^>^, ^ + A / <0, (32) 

and, similarly, for the third area 

v < ^, v + A / > 0. (33) 

The relations (30)- (33) are convenient that allow to easily estimate the size of epidemie 
even at the initial moment of time on the value of size v + XqIq- 

In summary, we shall note, that in view of complexity and an insufficient level of 
scrutiny of the mechanism of immunity, the corresponding mathematical models, and 
also, based on them, model of epidemies are, certainly, enough rough. Nevertheless, 
apparently, even the estimations received in the elementary models, in practice can be 
appear rather effective. 
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